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ABSTRACT 


The problem of Multi-user communication with Spread 
Spectrum Multiple Access over fiber-optic channel has been 
considered. Specifically Direct Sequence-SSMA situation is 
considered. Suboptimal receivers, namely, a correlator, and a 
cascade of inverse filter and a correlator are proposed for 
this problem in unlimited and limited bandwidth cases 
respectively. The latter case is introduced at the receiver 
(not over the channel) in order to reduce the thermal noise 
effects of the front end circuitry. The performance of the 
system is evaluated taking into consideration the effect of 
shot noise. The probability of bit error is evaluated using 
the Gaussian approximation f or the decision statistic* This 
approach allows us to see the performance dependence on various 
parameters, like optical power, APD gain, code length of PH 
Sequences# Bias .Resistance and number of users etc., more 
explicitly* The calculated results are presetted along with 
the raettits of simulation made for the unlimited bandwidth case. 
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CHAPTER 1 


INTRODUCTION ; 

1.1 THE MULTI-USER COMMUNICATION PROBLEM: 

The problem of two or more users communicating over a 
single channel has been considered in the past. Such a case 
exists in two ways. One is that a group of users is situated 
at one end of the channel, and all of them want to communicate ■ 
to a similar group of users at the other end of the channel j 

and vice versa. The second is that the users are geographically 
apart and are connected by a network of communication, each 
of them wanting to communicate ^"^any other. To achieve this 
‘Multiplexing’ is used in the first case and 'Multiple Access’ 
in the second case. A vast majority of the existing Multiuser i 
conmunica tion systems use ‘Frequency division Multiplexing* 

( FEM) or ‘Time Division Multiplexing* (TDM) for the first 
case and 'Frequency Division Multiple Access' (FEMA) or the 
'Tim© Division Multiple Access' (TIMA) for the second case. 

i 

[ 

Though known for quite some time the third method | 

'Spread Spectrum Multiplexing C SSM) or 'Multiple Access' (SSMA),! 

i 

is being considered as a viable alternative to the above ! 

■’ *" " * < ' ' ' ----- • • ; - I 

, ;■ ■ , .. 

methods during recent times. Unlike the other two methods j 

SS Modulation offers some additional advantages in addition j 

to providing multiplexing or multiple access capability. In 
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fact these advantages only motivate us to consider SSM(A) 
instead of TEM(a) or FEM(A). The advantages are narrowband 
interference rejection, tone frequency (ac or dc) rejection, 
antijamming capability, and antimultipath capability [l ], 

It also offers a certain amount of secrecy. Some of these, 
we show more analytically in Chapter 3. As a multiple access 
method SS modulation offers an important feature, namely that 

: 3 J 

of asynchronous multiplexing of the users’ signals. That is 
in a network users can be added or dropped with out affecting j 
the performance of the other users. Thus it obviates the need 

I 

of expensive loop synchronization as in the case of TEMA. 

Though this feature of asynchronous mode of operation is present 
in FEM& also, there the use of the channel bandwidth by a j 

particular user is permanately limited to a portion of it 

| 

while all users in SSMA simultaneously utilize the channel 

I 

fully to achieve the multiple access in addition to obtaining 
the above mentioned advantages [ 1 ]. 

I 

The technique of SSM(A) can be explained as follows. 

Each user is assigned a spectrum spreading sequence which is . j 
near orthogonal to that of any other user. Such spectrally 
spread signals are added and transmitted over a channel (in 
SSM). At the receiver this composite signal is correlated 
with each of the SS sequences thus recovering individual 
user signals. The actual signal is recovered while all other 
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interfering signals are spread out so as to cause minimum 
interference. The method is explained in detail in Chapter 3, 
The bandwidth of the spreading signal is very large when 
compared with the bandwidth of the information bearing signal 
it spreads. The higher the ratio of the two bandwidths the 
larger the advantage we gain by SS modulation. However, this 
is not the case with the Additive White Gaussian Noise (AV/GN) 
environment, where SS modulation does not offer any advantage 
over other schemes. Later we will have occassion to show 
this. 

The high bandwidth requirement of the SS modulation (or 
S£M A) if used in a data network can be met by optical fibers. 

In addition to this, many other advantages of the optical 
fibers [2 ] led people to consider it as a medium for different 
data networks [3 ],[4 ]• Both of them put together, SSMA in 
optical fiber communications seems to be a good alternative 
to the existing schemes. One such attempt is given in [3 ]. 
However, such a conclusion can only be made after analysing 
the performance of S3vIA in fiber optic communications. Unlike 
the traditional communication systems, where the AW3N is the 
fundamental limitation on ideal c omniunication optical commu- 
nication involves an additional fundamental source of noise, 
namely the shot noise. It is inherent to the photodetection 
process, at the receiver, by which the optical power is 
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converted to electric current. So a complete analysis is 
required to see the effect of such shot noise and the 
limitation it introduces into the system performance. This 
has not been considered in [ 3 ] . An attempt has been 
made here in this thesis to study this and arrive at the 
system performance of SSMA in fiber optic communication, 

1.2 OUTLINE OF THE THESIS: 

This thesis consists of five chapters. In Chapter 2 
we scan through the problem of communication over a fiber 
optic channel and the two approaches taken by people to solve 
the problem of detection at the receiver. 

In Chapter 3, the SS modulation and its capabilities 
are discussed. The multiple access situation and the 
performance limitation it introduces are also analytically 
presented. We also discuss briefly the different methods 
of SSMA. 

In Chapter 4, we propose the receiver and evaluate 
the performance of the receiver. Both the cases of unlimited 
bandwidth and limited bandwidth are considered. Inverse 
filtering is proposed as a possible solution to the problem 
created by the bandwidth limitation at the receiver. The 
system for unlimited bandwidth case is simulated on a 
digital computer. 
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In Chapter 5, we present the results of the performance 
predicted in Chapter 4, along with the simulation results. 

The chapter concludes with a critical evaluation of the 
proposed system and the limitations of it and with suggestions 


for further work 



CHAPTER 2 


DIGITAL FIBER OPTIC COMMUNICATION RECEIVERS 

This chapter begins with a description of a general 
model of the digital fiber optic communication system. The 
model suggests the different signal processing schemes that 
can be implemented. We then discuss briefly the wo rk done 

! • 

by others in this field previously. 

2.1 THE FIBEROPTIC COMMUNICATION SYSTEM MODEL : 

A general digital fiber optic communication system is 
as shown in Fig. 2.1. The information source puts out random 
binary data which are transmitted using a suitable pulse i 

I 

format. The typical pulse formats are PAM and PPM. They 
are as shown in Fig. 2.2. Other pulse formats can also be 
used. These pulses modulated by the data bits are fed to a 

1 

LED or Laser which illuminates the fiber accordingly. If the ; 
intensity of the light incident is directly proportional to j 

the signal pulse, it is called Intensity Modulation. This j 

2 

intensity modulated optical signal, can possibly be attenuated j 

- >■ «... . • V. - I 

and distorted while travelling on the fiber, before reaching ! 

I 

the receiver end, | 

t 

At the receiver this optical signal can be processed in j 

i 

two ways. One is to take the electromagnetic considerations j 

i 

j 

and process the light accordingly [6], Such a method calls 
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for coherently detecting the light. But these methods are 
complex and highly sensitive to allignment errors etc. The 
other method is to photo-detect the optical signal, i.e. 
converting the optical power into electric signal, prior to any 
signal processing. Such a scheme is known as ‘Direct 
detection’ method, A photo-detector, either a PIN diode or 
an Avalanche photo diode (APD) is normally used for this purpose, 
The former simply converts the optical signal into electric 
current without any amplification resulting in a low signal to 
noise ratio (SNR), The latter provides signal amplification 
resulting in a high SNR. In both cases they detect the signal : 
as a noise process whose average value is related to the signal. 
This noise is known as shot noise and is of fundamental 
importance to the fiber optic communication problem. What 
distinguishes it from the regular thermal noise is its dependence; 
on the signal. The average value of this process will be : 

t 

linearly related to the signal transmitted provided we assume 

I 

that optical pulses add linearly. Such a case has been 
considered by Person ick [5] and he shows the assumption to be 
true for multimode fibers. The photodetector has associated 
biasing circuitry and subsequent amplification of the signal 
also takes place. This front-end circuitry in addition to, 
possibly, distorting the detected signal, adds the thermal 
noise. At the circuit level the various noise sources can 
be as shown in Fig. 2.3, 
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The thermal noise source associated with is 
modelled as a shunt current source i^( t) . The amplifier 
noise source, referred to the input are i^( t) and e^( t) . 

They are followed by the noiseless ideal amplifier with gain 
A. The current source i ( t) 'the signal source* accounts 

o 

for the shot noise current. It is actually a sequence of 
randomly located impulses whose emission times are governed 
by the intensity of the incident light, >«,( t) .Specifically, 

k(o, t) 

i (t) =2 eg- (t-t.) (2.1) 

5 i=l 


Where ’ e 1 is the electron charge and g^ are the random gains 

i 

provided by the APD (i.e. corresponding to each primary 
electron a bunch of g^ secondary electrons is released). The j 
density function of g^ is usually complicated [15], but its | 

moments can be calculated knowing the ionization constant of j 

the photo detecting device. In case of a PIN diode 1. I 

t^’ s are the random instants of time when an electron-hole 
pair is generated due to the absorption of a photon. k(0,t) is : 

known as the count variable and is' equal to, .the number of 1 

j 

electrons emitted during the interval (0,t). Clearly k(0,t) is j 

j 

a random variable h It has Poisson statistics with density [6] | 


* A* 

p(k=n) = 


-A 


t 

where A= / A(t)dt 
0 


( 2 . 2 ) 

(2.3) 
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which means that the probability of emission of n electrons 
in the interval (o,t) is given by the above expression. 

We know that intensity X(t) of the optical field depends 
on the time varying optical power as 

A(t) = ~hy P(t) +A 0 (2.4) 

where n : photon counter efficiency 
hy : energy of the photon 
A- : dark current intensity 
and p ( t) ;is given by 

p(t) = Z d. tv(t-iT) (2.5) 

i= — o° H * 

where d^ are data bits £o,l £and h ( t) is the optical channel 
output pulse as shown in Fig. 2.1. 

T - data bit duration. 

If v Q (t) is the output of the amplifier of gain A, then 

k(o, t) 

v Q (t) = A 1 e g ± h fe ( t-t i )+A(i b ( t) + 

i A (t))* h fe (t) + A e A (t) (2.6) 

where h^ c (t) is the impulse response of the front end RC 
circuit, i.e. we have a composite signal which is the sum of 
a shot noise term, whose average value is the signal desired 
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(we will show "this latter) and three independent Gaussian nois€ 
terms. Now the problem is how should v ( t) be further 
processed so that the desired signal can be detected with 
reasonably good fidelity. Infact, certain amount of noise 
suppression can be done by properly choosing the front end 
circuitry itself. Vie will see the consequences of such a 
step later. Before that we go through the work previously 
done to solve the above problem. ; 

In the literature basically two approaches are taken to 
process v Q (t). One is to have v Q (t) optimally processed. For 
this the likelihood function is evaluated and the resulting 
structure implemented. Such receivers are obtained by 
Bar- David [7], Gagliardi and Karp [6], and Foschini et.al. [8]. 
Bar-David has considered the problem with no ISI and AWGN 
free situation. Gagliardi and Karp considered the problem 
with no ISI and additive Poisson Noise. The structure arrived : 

at by the former is a cascade of electron emission times i 

i 

estimator and a weighted photon counter. The Maximum ft | 

^Posteriori (MAP) receiver derived by Gagliardi and Karp finds a | 

| 

pair of weighted sums of the observed count components j 

and compares the sums to get the larger. j 

Foschini et.al. considered pulse dispersion also and j 

which is a cascade of tw> viterbi 


they obtained a receiver 
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algorithm (VA) receivers, one to estimate the electron 
emission times and the other to estimate the data sequence 
in the presence of ISI [9]. The impracticality of VA 
structure at high data rates led to the consideration of the 
possible use of Decision Feedback Equalizer (DFE) in its 
place. Such a subop timal structure vas obtained by Govind 
Sharma [10], However this structure also requires high 
sampling rates at the receiver which again limits the data 
rates. It also requires a high signal to noise ratio. The 
ISI problem is also considered by Dogliotti et.al. [ll] in 
which either a transversal filter or a Kalman filter is 
offered as the receiver structure. Shot noise was not 
considered in this, Messerschmitt [12] also arrived at the 
cascade of a forward filter and a DFE which minimized the 
mean square error. In the last two works performance 
evaluation in terms of probability of error was not done. 

The other approach is to subop tirnally process the v Q (t) 
by assuming a receiver structure which does not consider the 
discrete nature of the shot noise to find the optimal 

structure like the first approach. However, the complexity 
of the optimal receivers led to this approach. Thus the 
work of Personick [13] proceeds this way where he obtains a 
linear solution to this problem. He takes v 0 (t) and passes 
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it through an equalizer which eliminates the ISI comple- 
tely (zero forcing) at the sampling instants by selecting 
the output pulse shape to be a raised cosine* Unlike the 
above methods where the exact method (of finding the 
density of the decision statistic) is used to find the 
performance, here only the Gaussian approximation (or the 
signal to noise ratio approx.) is used to find the perfor- 
mance, i.e. the probability of error. He argues that for 
practical systems of interest results obtained from such 
calculation differ by about 2 db from the ones found by the 
exact method. 

In this thesis we attempt to solve our problem on the 
lines of second approach. The explicit statement of the 
problem is: How do we process v Q (t) and obtain the performance 
of such a processor when v ( t) is owing to the sum of K 
asynchronous SS modulated signals. The system can be shown 
as in Fig. 2.4. The suboptimal solution is sought for the 
reasons that in the case of AWGN its performance is reasonably 
good and the complexity of the receiver is also small. A 
fiber optic SS modulation case has been considered earlier in 
[14] with no multiple access and no ISI. The receiver given 
there does not suit for our case for reasons explained in the 
next chapter. Before going on to propose the actual receiver 
structure we discuss in the next chapter the various capabi- 
lities of SS modulation and SSMA techniques. 















CHAPTER 3 


SPREAD SPECTRUM MULTIPLE ACCESS COMMUNICATION SYSTEMS 

In this chapter we first present the various 
capabilities of the SS modulation. After that we give the 
reasons for choosing this Multiple Access method instead of 
TDMA or FDMA. Lastly we give the Multiple access system model ; 
(without any reference to fiber optic communication) and presen' 
the analysis for quantifying the multiple access interference. 

3.1 SPREAD SPECTRUM COMMUNICATION USING A RANDOM SIGNAL MODEL : 

The concept of Spread Spectrum (SS) communication and 
its Multiple Access method can be easily understood by assuming 
that the signals used for spreading of the incoming data 
sequences are random processes, which are not suited for system 
implementations. For this reason, in what follows, we will 
not give a too mathematical treatment. A more precise analysisj 
will be presented after this for a specific system model which ' 
is implementable* I 

Let c(t) be a zero mean stationary Gaussian random | 

process which has autocorrelation function R (t) given by | 

V 

R c (t) = (T c r X (T c - lx/ ), hi ^ T c 

T c 


= 0 


(3.1)' 
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The actual shape of R (■£) is not important as far as it is 

w 

approximately zero out side (o,T ). Its spectral density 
S c (f) is given by 

S c (f) = T c [ sinc(fT c )] 2 (3.2) 

Let s( t) be a binary data signal, of duration T and 

constant amplitude +A (A>o), We assume T»T . Consequently 

the bandwidth of c( t) is large compared to that of s(t)« 

Assume a noiseless baseband channel over which we send the 

SS modulated signal, s( t) c(t), and use a correlator in the 

receiver to recover s(t). The decision statistic is 

T T T 

2 = / [s(t)c( t)]c(t)dt = / s ( t) <?( t) dt = + A / c 2 (t) dt 
o o ~ o 

s( t) £ £a,-A 

Hrgodicity implies that with high probability 

r 1 /c T 2(t)dt = E[<At)] = R (o) = 1 
0 C 

Z = + AT 

Thus s(t) is recovered exactly in the absence of channel noise. 
However it is to be noticed that, for this recovery, the 
knowledge of c( t) at the receiver is necessary which : makes 
this conceptual model non implemen table. Usually c(t) is 



approximated by a pseudo noise (PN) waveform derived from a 
PN-sequence. We now consider the effect of interfering signal 
on this SS communication system which we claimed in the first 
chapter to be reduced due to SS modulation. 

1. Anti- Interference capability: 

Let the received signal be 

f ( t) = s( t)c( t) + A’ 
where A* - d.c* signal 

Now, ^ T 

Z = / s(t) c(t)dt + A* / c(t)dt 
0 o 

= + AT + A' T { r 1 / c( t) dt l 

T 

The interference is negligible if ^ / c(t)dtAf 0 
which is true if A’ as A since c( t) has zero mean and T is large 
compared to T * Thus it can reject the d.c. interfering 
signal. By using similar argument one can show that it has 
a.c. interference rejection capability by making use of this 
fact that 

i T 

~ / c(t) cos u Q t dt»0 
o 

where is the frequency of the a.c. tone* 

2. Anti-Multipath Capability: 

In a simplified base band model of a single multipath 
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channel the received signal might be 


r(t) = s(t)c(t) - f3s(t~ r)c(t-t ) , O < p < 1 
If s(t) = A for (o,T) and * > T 

C 

T T 

Z = / r(t)c(t)dt = AT - / c( t) c(t- ) dt 

o o 

«*1 - pR c (x)AT AT 

Thus if the multipath delays are greater ’than T, they 

v 

do not affect the performance of the system. 

3. Multiple Access Capability: 


Let there be K transmitted signals, each with a 
different process c k (t), 1 < k < K, The processes are 
mutually independent. The received signal is 

K 

r(t) = 2 s,(t) 

k=l K R 

The output of the ith correlator is 
T 

Z . = / r(t) c.( t) dt (3.3) 

o 

T K T 

=» / s.( t)c.( t)dt + 2 / s. (t)c k ( t)c. ( t)dt 

o 1 1 k=l o k K 1 

K 

= + AT + 2 + A / c ( t)c. ( t)dt 

i=l o R 
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Since 

c k (t), c i ( t) are independent zero mean processes 
r" 1 /^c^C t)c^( t) dt as 0 

Z j _ = £ AT. This demonstrates the multiple Access 
capability of SS modulation. 

In all the above cases, the key feature is that the 
c( t) has low correlation with the interfering signals. Thus 
It produces only a negligible change in the output of the 
correlation receiver. The integrator tends to 'smooth* the 
interfering signals, but it produces a large peak in response 
to the desired signal. Such an explanation can also be 
given in terms of frequency domain. 

In addition to the above it can be shown that SS 
modulation can also combat narrowband, partial band 
interferences [l}„ 

3.2 SYSTEMS OF SS COMMUNICATIONS: 

There are mainly three methods and some combinations of 
them to implement the normal SS communication system they 
are: 

1. Pseudo Noise Sequence modulation 

2. Frequency Hopping modulation 

3. Time Hopping modulation 

In the first method the spreading of the spectrum is achieved 
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by multiplying the data signal by a (usually) binary PN - 
sequence. This is also known as direct sequence SS modulation. 
This system is easier to implement. The second method uses 
a pseudo random assignment of the carrier frequency of the 
base band data modulated signal such that the signal in 
frequency domain sweeps the frequency band by hopping in a 
pseudo random way. These are relatively difficult to imple- 
ment due to 'the difficulty in implementing the frequency 
synthesizers with such a wide range. The third method, namely 
the Time hopping, transmits the signal only at selected 
intervals (varying) of time chosen in a pseudo random way. It 
has an advantage that the average power requirement is less. 
Normally a combinations of any two of the above seems most 
useful. Many such systems are given in [2.6] , [27], [28], 

In the next section we consider the first method, i.e, PN - 
sequence modulation with reference to Multiple Access 
situa tion. 

3.3 DIRECT SEQUENCE SPREAD SPECTRUM MULTIPLE ACCESS 
COMMUNICATION : 

Before going on to describe the SSMA, we give a 
heuristic argument for choosing SSMA as opposed to TDMA or 
FDMA. In an ideal condition all the three offer equal number 
of users for a given system bandwidth. However, when other 
aspects such as fading due to multipath, interference from 
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other communication sys terns, intentional jamming etc* affect 
the performance, the use of SSMA over the others proves 
useful. The other important advantage of SSMA over TDMA is 
that it can be operated in an asynchronous manner (i.e. 
users are asynchronous), a fact which makes us prefer SSMA 
to TDMA because of the costlier synchronizing equipment needed 
in case of the latter in a Multiple Access environment. 
Especially in data networks this asynchronous mode of operation 
can also be exploited in terms of reducing the complexity 

of network protocols. In fiber optic communication, the 
effects of nonlinearities of the transmitting and receiving 
devices reduce the performance of FDMA where as SSMA has the 
capability to reduce their effect [3], For these reasons the 
SSMA seems to be a viable alternative to the other two. 

However it should be noted that the advantages of SS 
communication are not present in some situations. For example 
in a purely AV/GN environment (no other source of interference ' 
is present) the use of SS modulation offers no advantage. 

I 

To see this, let the received signal be 

I 

1 

| 

r(t) = d( t)a( t)+a)( t) ; d( t) - message signal, + 1, 

N 9 

a(t) - PN sequence, w( t)-A’j®I zeromean variance =-£■ j 
+ 1 

As usual we use the correlator to recover the data ! 
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T 

Z = / r(t)a(t)dt = d T + w, (3.4) 

o 01 

E [Z/d Q = +1] = T i 

T T 

Var [Z 2 ] = f f a(t,)a(t )E[w(t*)io(t 1 )]<\t 1 dt 2 
oo ^ ^ • 

NTT N 

= 2 a / / a (t x )a(t 2 ) at 2 = 2^ T 

2 

P e = erfc*(d) where d 2 = = §- T 

Thus it has nothing to offer as an advantage over BPSK. modu- 
lation scheme in this situation. 

The effect of shot noise on the SS communication will 
be seen in the next chapter. We now give the multiple 
Access situation with its interference effects with no other 
additive noise sources. 

Binary Baseband SSMA Communication: 

In binary direct sequence form of SS modulation, a 
baseband signal of the form 

x(t) = X x-P T (t- jT ) (3.5) 

3 L c C 

is employed as the spectral spreading signal. In the above, 
sequence ( x j) is a periodic: sequence of elements ^+1, -1 j. 
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and P_ ( t) is a time limited signal [limited to (o,T )] 
l c ' c 

for which 

T c 

C 1 / Pt (t)dt = 1 (3.6) 

o c 

The normal choice for fibre optic communication would be 

P T (t) = 1 0 < t < T c 

c 

=s O otherwise 


P~ (t) is commonly referred to as chip waveform. The sequence 
A c 

x = (x^)» called the signature sequence when employed for a 
particular user signal, is periodic with period N* i.e, 

Xj s* Xj +N for each j. 

The binary data signal waveform d( t) is given by 


d( t) = z d, P T (t- ir) 


(3.7) 


£ = — 00 


where Pj(t) is non-zero ovdr [o,T] . d^ is the binary random 
sequence £+1, -jjcorr. to the user signal. When x(t) is 
used as a signature sequence for d(t), N is normally chosen 
to be equal to i.e. T is an integer multiple of T • So, 

for each data bit, a period of the sequence is modulated by 
each data bit d . Therefore the bandwidth of x(t) is of the 
order of N times that of d(t). With this, we now give the 


SSMA model. 
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As shown in Fig. 3.1, the SSMA has K asynchronous 
users dj(t) each with a random delay Each of these 

are modulated with synchronous (to the particular d(t)) 

, PN-sequences a^( t) and are linearly added to form a composite 
signal s( t) 


d,Ct) 


t 


a 

d^Ct -^> — * 


IP 


iVaJ!6 «* 


T-frlj 






(tecett 


*00 


Aa K Ct-rJ 


Fig. 3.1: Baseband SSMA communication system 
K 

s(t) = A Z d ± ( t- x i )a i (t- x ± ) ' (3.7) 

= . Z s.(t) 


The PN-sequences a^( t) are chosen so as to approximate the 
c(t) of section 3.1, so that their cross correlation and 
out of phase autocorrelations values will be small. At the 
receiver s( t) is ^ eC by noise n(t) (in traditional 

systems n(t) is AV/QM plus any other interferences present) 
which in our case will be sum of thermal and shot noises. 
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As shown in Fig. 3.1, the SSMA has K asynchronous 
users c^(t) each with a random delay Each of these 

are modulated with synchronous (to the particular d(t)) 
PN-sequences a^( t) and are linearly added to form a composite 
signal s( t) 
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Fig. 3.1: Baseband SSMA communication system 
K 

s(t) = A d ± (t- i i ) a i ( t- (3.7) 


= . Z s.(t) 
i 1 

The PN-sequences a^(t) are chosen so as to approximate the 

c(t) of section 3.1, so that their cross correlation and 

out of phase autocorrelations values will be small. At the 

c orrupted 

receiver s(t) is £ • by noise n(t) (in traditional 

systems n( t) is AV/CN plus any other interferences present) 
which in our case will be sum of thermal and shot noises. 
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For the moment, we neglect their effect and consider only 

the multiple access interference. It is to be noted that all 

the users are assumed to have equal amplitude A. In such a 

case the detection of chips of each sequence and thereby the 

the data bits is not possible as the signal to interference rati 

is far less than unity • For this reason the receiver 

proposed for the SS modulated signal in fiber optic communica tio: 

in [14] is not applicable here. The equal power situation can 

be generalised to different amplitude levels (more 

appropriate for Multiple Access situation v\here the 'near-far* 

problem comes into picture) which modifies s(t) as s(t) = 

K 

E A.d.(t— t . )a . ( t— t). It is also assumed that tls are 
i_l i 1 1 1 

uniformly distributed over (o,T), For analysis purpose, we 
consider the receiver to be . synchronised to user one signal 
that Tj, can be assumed to be zero as shown, 

K 

r(t) = Ad 1 (t)a 1 (t) + A L d ± ( t- x ± )a ± (t- 

We assume the receiver to be a correlation (or matched filter) 

receiver. The output of this receiver is 

T T 

= / r(t)a (t)dt = A / d x (t)a 2 ( t ) dt 
1 o 1 o 

K T 

+ A E / d.(i> T i )a.( t- T i )a^(t)dt 
i -2 o 

= Ad ljQ T + A 1 (d^ S x a i (t- x i )a 1 (t)dt + d i>Q 

T 

J * 1>* X j^) ( "t) dt) 


( 3 . 8 ) 
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which is clear from the figure shown in Fig. 3.2. The 
integrals in eq. (3.8) are given by 

R i,l = /a j_( t- T i )a 1 (t)dt 


T 

R i,l = f a i^ t ~ T i) a x^ t ) dt 



Figure 3.2: AAsynchronous data sequences. 


which are known as continuous time partial cross correlation 

I 

functions. These we re introduced in [16] and used in the 
analysis of binary direct sequence SSMA communications [l6],[l7] 
Their effect on the receiver output is via the quantity. ; 


^k,l ^ d i* x ^ ~ ^ 


.-1 


^ d i,-l R i, 1^ 


+ d 




R i,l 


(r) ] 


(3.9) | 

/ 

[ 

which is referred to as the (normalized) multiple access interf^ 
ence at the output of the first receiver due to the ith signal. | 
The vector d- is the vector (d. , , d. ) of two consecutive | 
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data bits of the ith data source. Therefore can be written 
as 


A 


l 


AT 



K 

Z 

i=2 


fi^Ui^i)] 


(3.10) 


It is clear for K=1 (i.e, no Multiple Access), the summation 
term is not present. In that case, in the presence of 
AWGN, t ^ is the sufficient statistic [18] and we decide the 
first user data bits by comparing with zero for optimum 
detection. For K>1, it is not the sufficient statistic and 
a more complex receiver can provide a statistic which gives 
near optimal decision. However since correlator receivers are 
relatively single to implement, the vast majority of direct 
sequence SSMA systems employ correlation receivers, even 
though they are subop timai. This is because of the relatively 
small performance improvement that can be obtained with a 
more complex receiver. 


The usual performance criterion, the probability of error 
can be evaluated using two criteria. One is to find the 
maximum value of the multiple access interference (worst case 
value) and find the P in the presence of AW®. Such a method 

V 

is taken in [16] and expressions for worst case value of 
multiple access interference (MAI) obtained. The other is to 
consider the mean square value of MAI which leads to 
signal to noise ratio as a measure of performance. Vi/e now 
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proceed to obtain its value. The reasons for taking this step 
will be clear later. In general the first one is not preferred 
because the worst case MAI occurs very infrequently. 


Let, 

d = ^ d 2,-l ,d 2,o ,d 3,-l ,d 3,o*" d k,-l ,d k,o^ 
i = ^3* • • ,x k^ 


we can assume d, t to be independent as they arise from 
independent physical phynomena. We also assume d^ ^ and d m n to 
be independent for all i^m and j^n, and so are the x-, x- 
for i^ j . From these independence assumptions, we see 


Etl i,l (d i, x i^ = °* Since P(d i,-1 = +1) = p ( d i, 0 ) = 1/2 ’ 
the variance of I. , (d. , x.) is given by 

1 y X X X 

4 ,x = v « Ch.iK’ = r 2 / r 1 [Sf > 1 (t)+ «? tl (»)]dt 


Defining the quantities 
we have 



(3.11) 



4,i = r3 < m i,i + “i.d (3 ' 12) 

The quantities m. ,, m. , are analysed and closed 

A 1 j X X j X 

form expressions are obtained in [19]. Using those expressions, 



can be written as 
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0 i,l ~ X t 2 ^i,l^ + 


where the parameter p. . (n) is given by 

J 


N-l 

v>* An) = s c A{)CAi+ n) 

1,3 {=1-N 1 3 


(3.13) 


(3.14) 


where C. H), CAl) are the aperiodic autocorrelation 
' J 

functions [16] of the i and j th signature sequences respectively. 
They are given by 
N-l-/ 


C,(/) =2 a, 


U) 


3=o 


a . 


(i) f 

j+Z , 0 < Z<N-1 


= 2 1-N <-i< 0 

j=o 3 ~< J » 


= 0, K| > N 
More generally. 


V * a oo a u ) 

j= 0 j J - 1 - 4 ’ 

0 < £ < N-l 


V +/ a< k > a U) 

3=o J J 

O 

V 

vl 

H 


0 , U| >N 


(3.15) 


and C^Z) = 
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where is the jth chip element of the ith sequence. Note 

J 

that C ± (l) = C ± (-£ ), l> 0. 

Thus from eqs. (3.13) to (3.15) it can be seen that a? • 

±9 J 

for all i and j (Mj) can be found knowing only IC Autocorre- 
lation functions C^^) of the PN sequences and K 2 cross corre- 
lation functions ^(X) are not required. With increasing 

0 

sequence length N, computation of erf 1 takes more and more 

1 j J- 

time. So some bounds are considered in [20] to minimize 
this time. 



9 

Once knowing the of . and the (possible) AWGN variance 

r, r 

we can find the signal to noise ratio as 


d 2 


a 2 t 2 


a 2 t 2 v J1 , 
A I £ O * -j + 

Mi 1,1 


rt 2 T 
0 th T 



a2s 4,i + 4/ t 


and the approximate probability of error (suggested by pursley 
[16]) is P tt <=* erfc*(d). The accurate evaluation, of P Q can be 
found in [2l] and in that the author shows this signal to noise 
ratio approximation is very close to the accurate value when 
the number of users is large and the users are asynchronous, 
with high signal to noise ratio. We will# examine the accuracy 
of this approximation in Chapter 5 by performing simulations 
of the system. 


Before concluding this chapter, we mention briefly 
about the selection of the above PN- sequences. For a given N 
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of reasonable length the Max. length sequences are relatively 
small in number. For example for N=511 they are 24. For 
N = 1023 they are 30. So as an alternative to these, we 
consider other sequences with these periods. Two of them 
are Gold Sequences and Kasami Sequences [22]. For every 
N = 2 n -l, n /o mod 4 there exist 2 n +l Gold sequences with 
period N. In Kasami Sequences there exist 2 classes, one is 
the smaller class and the other is the larger class. In the 

smaller class there are 2 n ^ 2 sequences with period N = 2 n ~l. 

2 n ' 2 ( 2 n +l) 

In the larger class there exis ^sequences with period same 
as N, Earlier, the»j sequences were found for the following 
reason.. When synchronous users are present, the Multiple 
Access interference (not the quantity defined above) depends 
on the periodic cross correlation of the user PN-sequences [22]. 
The Gold and Kasami sequences were obtained in an effort 
to reach the lower bounds of their cross-correia tions in 
such a situa tion .With M. sequences such bounds are mot for 
only a few selected pairs of sequences known as preferred 
pairs. Quantitatively this bound [22] is given by 


.(>0 > 2^ n+1 ^ 2 + 1 , 

J 

n odd 

2 (ih 2)/2 + lf 

n even 
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Fox Gold sequence the equality is satisfied provided 
n is ^ o mod 4. 

In the present asynchronous case, however, no such 

2 

sequences optimal with respect to <• interference variance 
are known till now except for a class of M.-L sequences of 
length 31 to the best of this author’s knowledge. So, for 
the problem at our hand, we have chosen Gold sequences 
considering their number. The selection and generation of 
the Gold sequences axe given in Appendix A. 



CHAPTER 4 


FIBER-OPTIC SSMA RECEIVER 

In this chapter we propose the receiver for the problem 
stated earlier in Chapter 2 and analyse its performance for 
various cases. In the first section, the case of unlimited 
bandwidth will be analysed. The next one is to assume an 
equivalent band limited channel and see how the proposed 
receiver for the previous case be modified to get satisfactory 
performance. This will be discussed In the second section. 

4.1 SSMA IN UNLIMITED BANDWIDTH CASE: 

The system model is redrawn with the proposed receiver 
as shown, for this case in Figure 4,1. The users* data 
sequences are d^(t), d^t)... d fc ( t) . For analysis purposes, 
only the detection of user 1 sequence, namely d^( t) will be 
considered as all the others will also go through same 
receiver operation. It is assumed that perfect synchronization 
exists between each pair of $S modulator and correlator 
corresponding to each user PN-sequence. The validity and 
consequences of such an assumption are discussed in Chapter 5. 
Since we are considering user 1, we make =* o without 
loss of generality. Each of 2 < i < K is uniformly 

distributed in (o,T)* where T is the data bit duration. The 
user sequences are SS modulated by the respective SS carriers 
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a 2 ^ ^ * • a j c ( ancJ summed up. This composite signal 

is fed to an LED (or Laser) which drives the fiber with 
an intensity modulated optical signal. The optical power of 
that signal is linearly related to the input driving current 
of the LED over a wide range of input current. 

In the fiber of unlimited bandwidth, the optical 
signal gets only attenuated while travelling, before arriving 
at the receiving end. At the receiver, prior to any signal 
processing, we let the optical signal be directly detected 
using the photodetector, either an APD or a PIN diode, as the 
case may be. With the former there is an amplification of the 
signal, though with an increased noise, in the process of 
detection whereas in the latter such an amplification is not 
present. Thus the signal to noise ratio (SNR) is larger 
in the former case than in the latter. The optimal 

amplification required can be derived. The detected signal is 
amplified and correlated to get signal statistic and is 
compared with threshold to decide the data bits. To see the 
operations, which the signal undergoes, we consider the 
equivalent models of the detector and amplifier and the 
correlator. They are as shown in Figure 4,2. 

The photodiode is modelled as a current source i ( t) 
which accounts for the shot noise due to the emitted inpulses 
of electrons at random time instants. represents the 




M.OPEL OF DETECTOR , AMPLIFIER AMP COfc,&£kAro&, 
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dio Co junction capacitance* is the bias resistance with 
its associated thermal noise current source i^( t) with 
power spectral density (p.s.d,) Am^Hz. R^ is the arapli- 
fier input resistance with the associated noise source i^( t) 
with p.s.d. Amj?/Hz. is the Amplifier input capacitance 

e^( t) is the series noise source with p.s.d. Sg volt^Hz. Other 
models of amplifier can be taken, but the basic approach is 
the same. 

The shot noise current source i g (t) is given by 


k(o f t) 

i s ( t) = E e g i 6 (t~t ± ) 


(4.1) 


i=l 

where g^ - Random gain of the APD 

t^ - random time instants at which electrons are 
emitted 

k(o,t) - count variable representing the total number of 
electrons emitted in )(o,t) 


v„(t) = i 5 (t) * h f .(t) 


f.' 


k ( o 1 1) 

£ e g* h f (t-tp 
1=1 


(4.2) 


where h^.(t) = impulse response of the biasing plus amplifier 

circuitries with time constant RC where 
R = r/r b ; C = C A + C d , 

i.e. h fc (t) = ^ e _t/I1C 
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The means and autocorrelations of i (t) and v ( t) are 

V 3 

given by [6] 


E [i s (t) ] = .1 e <g>X(t) + A e<g> (4.3) 

s hr ° 

where, 

t) - quantum efficiency 
h - planfck's constant 
y - frequency of the light used 
e - electron charge 
<g> - Average value of the APD gain 
A q - Dark current intensity 

A(t) - is the time varying intensity of the incident optical 
field which depends on the composite SS signal. It is 
given by 

A(t) a $ Jj(t) a . -if (K + d i (t-T i )a i (t~T i )) 

(4.4) 

where p(t) is the received optical power, b^ is the' bit .energy 
parameter which is related to the signal power at the trans- 
mitter as given below. T is the data bit duration. 


At the transmitter we have the composite SS signal 


a s 


s 



d i< 
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It is assumed that d i e|+l,-l^ , c £+1,-1 1 and 
T o 

/ ar.(t)dt « T. 

o x 

To this we add a d.c. signal of K units to get s( t) = K+s^ n (t) 
which varies in [0,2K] unlike s^(t) vhich varies in [-K,K]. 

We use this signal to 5 drive the LED/Laser. Since the input 
signal to an LED is in terms of its drive current > we 
assume 

K 

I(t),= l[K + E d..(t~ T i) a i( 1 >- \)]> 

I is the current in Amps. The emitted optical power is related 
to this current by 

p(t) » Q I(t) where Q is the responsivity (watts/Amps) 
of the LED/Laser 

it 

a QI [K + £ d. ( t- t J a 4(t~T^)] 

i=l 1 1 1 

So letting b £ « QIT (we multiplied by T since bg is energy 
parameter) , we have 
b F 

p(t) = (K + E d^t- T j.) a i( 1> \)) 

T is retained in the denominator only to give an idea of the 

bn 

power required* since -y has units of power. We now 
continue to find the expression for the mean of the shot noise 
process v s (t). 


(4.5) 
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E [v s (t)j =e^g)A( t) * h fe (t) 
neglecting the dark current intensity. 

Ri( w ” E [‘A)i,(V) = « E[g?] X(tj,)s 


+ E[i s (t x )] Efi s (t 2 )] (4.6) 

Let E [g|] = <g 2 > 

E [g] =f <g> 

H v ( ^ , t 2 ) = E[v s ( t 1 ) v s ( t 2 ) ] = e <g 2 > f\( t) .. 

h fe ( t l ~ T ) ~ T ) dT + ECVsC^)]* 

E[v s (t 2 )] (4.7) 

Thus it is seen due to the non-sta tionarity of the process, 
the mean and autocorrelation functions of i ( t) and v ( t) 
axe time dependent. We need these expressions in finding 
the SNR at the output of the correlator. In what 

follows, we consider that the fiber is dispersionless in 
time (i.e. causes No ISI) and analyse the system for two 
cases. The first is that the front end circuitry also does 
not spread the pulse. R and C are adjusted so as to make 
this condition prevail. In the second case, we consider 
the pulse spreading due to the front end. The reason for 
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introducing such deliberate dispersion will be clear when we 
analyse the first case. We now go ahead to do this. 

Let v Q (t) be the output of the amplifier with gain A, 

It is given by 

v Q (t) = A v s (t)+/y g (t) + Ae A (t) (4.8) 

where v g ( t) = [i b (t) + i A ( t) ] * h fe (t) (4.9) 

is the Gaussian noise contribution. It may not be white 

depending on h^ e ( t) 
or 

v D (t) = A i s ( t )* h f C ( *)+*( V t)+i A^ ^bf «Ct)+ Ae A (t) 

we now assume h^ ( t) is impulsive and is given by h^ e (t) =R6(t). 
Therefore 

v Q (t) = AH i $ (t)+AR(i b (t)+i A (t)+A e A ( t) 

We can write i (t) as sum of an average component <i (t)> 

S 5) 

and a noise component i ( t) (to minimize notation we retain 

o 

i s (t)) with zero mean. The autocorrelation of i s ( t) will 
therefore be 

Rf( ti» t 2 )- e <9 > X ( t x ) 6 (4.10) 

<i s (t)> is given by <i Q (t)> = e<g>p(t) + e<g> X Q 

The spectral demities of i b (t), i A (t), and e A (t), given 
by S ,S A ,S £ are 
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e _ . C _ 2k© . q _ 21c© 

S l> - R b ' S A - ' S E = — 


m 


where, 

k = Boltzman constant, Joules/*K 
© = Temperature in 

g m s* Transconductance of the amplifier 


mhos 


V/e assume the amplifier gain to be high enough to make any 
subsequent noise additions negligible. Let i be the statistic 
we obtain by sampling the correlator output at the end of each 
baud T. 


• • £ » / v o (t)a 1 (t)dt 


= A / [(<i s ( t)>+i b (t)+i A (t)+i s (t))Rfe A ( t)]a 1 (t)dt 

(4.11) 


( X* + Ah + ^ + ^ ns + ) 


w 


where 


jL = AR / <i (t)>a.(t)dt - ‘signal’ component 

s 0 

# T 

A b y a = AR / i b ^ a ( t)a^( t)dt - ‘Gaussian Noise’ components* 
° In this case they also are 

white. 

* e = AR / i ( t)a ( t)dt - Shot noise component 

n s . 5 jk 


* A / e A ( t)a v ( t)dt 


’White Noise’ component 
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We now find the square of the mean and variance of 
take their ratio to get the 'Signal to Noise Ratio’ i.e. we 
are assuming the statistic to be Gaussian random variable. Such 
an assumption makes the performance calculation simple and allows 
us to see how the performance depends on various parameters [15], 
In ordinary binary data transmission over fiber optic links 
with added Gaussian noise at receiver, such an assumption 
predicted the signal power required reasonably, close to the 
value found by exact calculations [15], [23], We rely on this 
and go ahead and see whether any of the additional features of 
S3MA make this, assumption more valid. In addition we assume ; 

that both shot noise and thermal noises are independent of 
each other and also the shot noise to be independent of the signal 
component <i (t)>. Such an assumption is normally made with 
the SNR calculations. Therefore, 

T 

i) A s = AR / (-J- e<g>p( t)+e<g> A o )a x ( t)dt 

„ T b P K 

= AR e<g> (K+ Z d.( T.)a ( t- t )) 

n » o i=l 1 x 

T 

a 1 (t)dt+ AR e<g>A Q / a 1 (t)dt (4.12) j 

’ | 

we have i 

T ! 

/ a. ( t)dt = C/N.T 
o x 



where 


C = average of value of the (sum of all binary chips) 

PN sequence. For ^sequences C = 1, and for Gold 

sequence, C = or -1 or, -(2^ n+1 ^ 2 +i) 

where the sequence length N = 2 n -l, n odd. 

l s = AR e<g> e . ~| • AR e <g> ^ 

/ T [ ai ( t)d x ( t)+ 2 d i (i>- T' i )a i (i> ^)] a i (t)dt 

o i=2 ^ 

+ e<g> ARX q . CT/N 

Usually A is 10 r /> of the signal intensity. Since now it is 
divided by N its effect will be negligible 

b E 

X s = AR ^ e<g>b E . KC/N + ^ AR e<g> ^ 


X X 4 

[ ./ d.(t) a 2 ( t )dt + / 3l (t) ■ 

o x o * 

f 

^2 d-^t- T i ) a i ( t- dt] 

l K ( 

= A f b^CC/N + A» b E [d x ^ 0 + y . E 2 \ d 

>1,-1 / i a 1 (t.T i )a 1 (t)dt + 

T 7 

o 

d i,o ^ a i (tr ” 3 

(4.13) 


/a 2 (t)dt = T . 


A ’ = AR s <g) 


Since 
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The first and the second terms in brackets indicate the signal 
and interference from . other users respectively; 

As far as the bit energy bg is concerned, both signal and 
interference are weighted by the same value. So increasing 
bg does not increase performance from Multiple Access inter- 
ference point of view. We take care of the first term in 
eq. (4.13) in the threshold setting (as it is a constant term for 
each correlator). The 2nd and third term appear as again | 

signal and noise terms i.e. 


i ’ e * *s “ A ’ b E d l,o + A ' b E’ T l d i,-l R i,l^ T i^ +d i,o R i,l^i^ 


= A 1 b c [d, + I I. (d., r. 

£ JL f Q 1 1 1 


)] 


(4.14) 


The interference term I.(d., x.) is familiar to us in Chapter 3 
where the mean and variance of it are given. The variance is 
given in eq. (3.13) 


E £ *si d l,o^ “ A ’ b E d l,o 
Since I i ( ^ i ) has zero mean and 

Vax [l s ] = (A') 2 b| J 2 o| (1 
ii)/ ns = AR / i s (t) ai (t)dt 


(4. 15) 


(4.16) 


E = 0 


by our zero mean assumption 
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T T 

Var [l nsl = a2r2 f S R i (t 1 ,t 2 )a 1 (t 1 )a 1 (t 2 )dt 1 dt 2 

= A 2 R 2 e 2 <g 2 > ^ 6(t r t 2 ), 

a l^ t 2^ dt l dt 2 by ec I*( 4 * 10 ) 

= A 2 R 2 e 2 <g 2 > /p( t 1 )a 2 ( 

f 0 

b- T K 

= A 2 R 2 e 2 <g 2 > ^ / (K+ S^d ± ( V i>T jL ))dt 

* ( 4 . 17 ) 

a 2 (t) = 1, te (o,T) 


Thus SS modulation does not have an additional advantage 
over other schemes in limiting the shot noise power* This is 
due to the fact that this shot noise is also white except for 
its non-stationarity and that SS modulation does not offer any 
advantage with white noise. Proceeding with the above 


b c T 


Var [J( ns ] =; A 2 R 2 e 2 <g 2 > JL (KT+d lj0 f a 1 (t)dt 

K ( T T y 

+ E 4 d. , / a.( fc- ‘tjdt + d. / a.( t-Odtv) 
i=2 L 1 »“ 1 o 1 1 1,0 i 1 1 J 


( 4 . 18 ) 
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Thus the variance of the shot noise term depends on all 
the other user data bits (d. ,,d. ). Our approach is to 

X f JL X 9 O 

consider its worst case value. Since for either M-sequences 
or Gold Sequences the maximum average value is always -ve 
(-1 for M-sequences) we can assume d^ d ^,o ^ or 
i^l. This makes 




A 2 R 2 e 2 <g 2 > $ £ (M + F + J 2 


- / a.( t- *c i )dt) 
o 


- A" b E (K + d lj0 < 4 - 19) 


where, 

A* ' = A 2 R 2 e 2 <g 2 > ^ , 

Ch - average value of ith sequence 

Since C^/N «1 we remove the explicit dependence of 

Var[Jt 1 on d, ( essentially same for d. = +1 or -l) 
n s J ipO m f v 

and assume it to be —1* In addition we assume G^C^C^C 
which is the worst case value in ^sequences or Gold sequences. 
Therefore, 

Var [A ns ] - A" b E (K + KC/N) 


(4.20) 
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Worst case C/N is 6^ for Gold sequences (and much less for bA~ 
sequences) and so the fractional increase (or change) of 
Var ^- s re l a "tively small due to the above approximations, 

iii) If We let = Jt a + = AR f (ifot t ) +i A^ t ^ a l^ t ) dt 


EC Ag] = 0 


Var [£ J = Var [jl 2 ] + Var ] 


[ i^(t), ij D ( t) are independent] 


A 2 R 2 / f T R a ( t lf t 2 )a 1 (t 1 )a 1 ( t 2 )dt x dt 2 


+ A 2 R 2 / / R b (t 1 ,t 2 )a 1 (t 1 )a 1 ( t 2 )dt x dt 2 


T T 


= A 2 R 2 [ / / S a 6 ( t i** t 2 ) a 1 ( t 1 ) a 1 ( t 2 ) dt 1 dt 2 


T T 


(4.21) 


o o 


I T 


+ f S s b 6 ( t 1 -t 2 )a 1 ( t 1 )a 1 ( t 2 )dt 1 dt 2 3 


o o 


A 2 R 2 T (S a + s b ) / a 2 (t)dt = T 

o 


= A 2 R 2 T (^P + ) - 2P . A 2 R 2 T (4.22) 

k A H b 


Here we observe that as and Rg increase variance of JL due 
to Gaussian noise decreases. This is where the deliberate 
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introduction of ISI is made at the receiver to reduce thfermai 
noise effect,' by increasing and R b (or equivalently R) . 
Before doing such a thing, we continue now with our analysis. 


T 

iv ) Jt w = A / e. ( t) a.(t)dt 
o 


= O 

T T 

Var L* w ] = A 2 / / S £ 6 ( VV a l ( t 1 )a 1 ( t 2 )dt 1 dt 2 


« 2 t 2k© 

i7“ 


(4.23) 


Thus by eqs. (4,16), (4.20), (4.22) and (4.23) 

Var DU = Var U s ] + Var [ i ns ] + Varfj 1+VarUJ 


= (A') 2 b|^ 2 aJ jl+ A-'b E (K + ^) 
a 2 r 2 t . ^ + a 2 t . 

R ' g m 

E 2 U] = E 2 [i s ] = (A«) 2 b 2 

2ron (A») 2 b| 


K 


KC 


(4.24) 

(4.25) 


SNR 


=> rJE-i a 

var[xj 


(A 


) 2 b| a f t i+ A ' ' b E (k+KC/N')+A 2 R 2 T. 


t 


+A 2 T. 

9 m 


(4.26) 



51 


K 

Z 

i=2 


2 

i,l 


1 

(K+KC/N)+ 

(A‘) 2 b E 


2kQ.A 2 R 2 T 
(A* ) 2 b 2 E 


(l/R+R 2 /g m 


Thus it is seen as b E ^bit energy’) increases shot noise 
and thermal noise effects reduce (latter’s more rapidly). 
Also as R increases the effect of e^t) also reduces as seen 
in the denominator. Beyond a certain value, any increase 
in bg does not improve the performance as the effect of 
interference is unaltered with bg. That is the performance 
‘levels off 1 after a certain bg, which depends only on the 
cross correlations of the individual PN- sequences, The 
greater the number of users, the larger is the shot noise. 

We observe later that this effect is small compared to the 
interference from other users directly. 

Substituting for A* , A' ’ 


SNR = -n> — 5 ^ — 

r of 1+ (ta <a% . fe-tK, |£) + (_^) 2 . fS? 

-_2 1 >- L M iN w 


e <g> 


(l/RHT 2 /g m ) 


<g 2 >=[k 1 <g>+(2- ^■)(l-k 1 )]<g> : 




V*- 


.2 


( 4 . 27 ) 


The use of APD with mean gain <g> is explicit now. <g z > and 
<g> are related by 




( 4 . 28 ) 
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where k^ is the ionization ratio of the photodiode. For 
silicon diodes k^is 0,01 to 0.1 (smaller the better) and for 
Ge diode it is approximately 0,5, The thermal noise 
decreases with increasing AP^ gain and hence the performance 
improves. However this is not mono tonic because increasing 
<g> increases the shot noise and so after a certain value 
of <g> the total noise increases with increasing <g> and 
performance falls down. By letting the total noise dependent 
on <g> equal to Nj, we can optimize <g> to have Nj take its 
minimum value. Therefore, 



(4,29) 


= ^i< s h( 2- IrW?) 2 


U/R + R- 2 /9 m ) 

b E 


letting G = <g> 


N t = A 1 (k ] ,Gf(2- ^)(l-k 1 ))+A 2 . 


where. 


A x - ^ • k « * > 


A 2 “ ( r p) 2 * 2 2^2 * + R 2 / g m^ 

G t) — 



53 

differentiating w.r.t. G and setting the result to zero 
we have 

r i 2A o 

d - A i • ~2*~ ]“ — -r = 0 

a 1 x or 

or 

A 1 Ki g3 +A 1 (R. 1 )&.2A 2 = 0 (4.30) 

The solution of the equation (4.30) provides us with the 
optimum gain <g> opt. 

In digital communication systems, the accepted 
performance criterion is probability of bit error. We 
obtain it simply due to our assumption of A as Gauss, random 
variable, as 

P e = erfc* (SN R 1 ^ 2 ) (4.31) 

where, 

erfc* (x) = f e~ t / 2 dt 

X 

We give the results of this criterion in the next 
chapter where we observe the variation of P e with b £ , K, 

<g>, N and R. The variations 'due to the last parameter, R 

motivates us to consider the case of ISI introduction by 

2 1 ^ ^ 

case in the 

next section. 
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4.2 BANDWIDTH LIMITATION AT THE RECEIVER: 

As we saw, the bandwidth limitation at . the receiver 
is caused by introducing a high impedance front end after 
the photodetector. Such a thing causes the data pulses to 
overlap and to introduce inter symbol interference. This 
problem, namely that of bandwidth limitation can even be 
generalized to the fiber channel also, in which case the 
pulses get overlapped in the optical format itself before 
arriving at the receiver. Such a condition is as shown 
in Figure 4.3. 



P<OC&3l$4* [ 


Fig. 4.3 

.The optical intensity X( t) at the input to APD is now 
modified as 


^(t) = hy "T ( K+ | d i (*- ^i)^h ch (t) (4.32) 

and, average value of the shot noise at the output of amp. is 

for, 

V 0 (t) = gy-T \))* h ch (t)« h f ,( t) 

• V hy “§ (K+ | d i (t>- 
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where h^(t) = h c ^(t)* h^ € (t) is the equivalent base band 
channel impulse' response. 

Now how to process the detected electric (composite) 
signal to decide on d^(t)? One approach is to take up the 
Gaussian noise method where in the signal is distored by 
additive white Gaussian noise and Inter-symbol interference [29], 
The interference from other users and the shot noise can be 
approximated as Additive and Gaussian, The other is to 
optimally process the detected signal with associated shot 
noise, AWGN and ISI [8 ], However, for reasons explained in 
Chapters 2 and 3, we tend to choose the suboptimum receiver. 

Thus we have taken up the first method where only a signal 
embedded in white noise is received after distored by ISI. 

So we proposed the matched filter equalizer structure for the 
problem at hand. It is as shown in Fig, 4.4, i,e. we use 
matched filter- dec is ion feedback equalizer structure to decide 
on d^(t) and obtain its performance. But such an approach 
immediately warns of the low signal to noise ratio obtained at 
the 0/P of the matched filter. The low SNR is due to the 
fact that all the SS modulated data sequences passed through 
the band limited channel by which they are severely distored 
and. the cross correlations of which are no more small but 
same as the order of the autocorrelations of themselves. 

In view of Chapter 3» this can be explained as follows. There we 
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have shown that the cross correlation of two sequences is 
dependent upon their auto correlations. The PN sequences, being 
Pseudo random have a very small autocorrelation outside the 
chip interval. Thus their cross correlations are also small. 

But now the PN sequences are so much spread that their auto- 
correlations will not be any more near- impulsive. Hence the 
cross correlations also will not be small. And in this low 
£>NR environment the performance of DFE is very poor. So such 
a structure cannot perform satisfactorily. For illustration 
we have shown in Table 1 autocorrelation and the correlation 
values of a PN sequence with another PN sequence with and 
without bandl imitation. The resulting low SNR is as expected 
above. This is infact the reason, historically, for using 
S3MA only on a channel with available bandwidth [24].- Otherwise, 
the spreading of the data or N has to be decreased accordingly. 
Infact, it is for the same low SNR reason, that chip level 
detection of sequences is not done in S3MA. 

So, to solve the problem at least partially, without 
losing the advantage we saw with high impedance front end. We 

took the following path. We assumed that the fiber does not 

* ’ / 

introduce any ISI. The possible 131 is only caused by the 
high imp* front end. For convenience we repeat the receiver 
block diagram in Fig. 4.5. 
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Fig. 4,5: SSMA in band limited case 

Since, we already chose the correlator for the SSMA 
case we retain it as it is. Now to make the SNR high at the 
output of the correlator we need to nullify the 

adverse effect of the h^ e ( t) but retain the advantage with 
high R. From the frequency domain point of view it is clear 
that h^(t) alternates the> high frequency components* So 
the siiqplest solution will be to follow it up with an 
inverse filter which accentuates the high frequencies, to 
make it explicit we seperate h^ e (t) into cascade of two units 
and put an inverse filter to the second one. It is as shown 
in Figure 4.6. 



Fig. 4.6: Inverse filter solution 
1 -i/RC 1 

The Fourier transform of ^ o ' =* RCjuh-1 


(4.33) 
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We choose F.T. of the Inverse filter impulse response 
to be (RCjoH-l) such that the cascade of the two filters 
produces an impulse response of R6(t). Physically the 
inverse filter is the direct addition of the output of h^Ct) 
with a differentiated and amplified version of itself# Thus 
at the output of the inverse filter (or at input to the 
correlator) the filtered versions of shot noise (including 
signal) and thermal noises i^( t) and i b ( t) are spectrally 
restored back to their original shapes. So the analysis of 
the unlimited bandwidth case is directly applicable. However, 
the white noise source e^( t) is affected by the inverse 
filter rather harmfully. In frequency domain, the e A (t) with 
a constant p.s.d. of S £ is now modified by as 


% B -%. |H lnv ( f )| 2 = SE CHW +l) =2|(RW +l ) 


fA-ZA) 


Its autocorrelation function is given by 

r (O . aa ( #(t) - r 2 c 2 -4«<t» < 4 - 35) 

G y m d 

If we let c (t) be the Gaussian noise corresponding to © A (t) 

and JL , the statistic contribution duo to it 
/ 00) 

/ ou = ;\ (t)a i‘ t)dt 


* 0 
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Var[ 
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T T 

SI 

o o 


R e (Y) a 1 (t 1 ) a 1 (t 2 )d t 


l dt 2 


II TT 

/ I S( VV a l (t l )a l (t 2 )dt l dt 2 + I I 


o o 


0 o 


2 2 ^ — ^) ) 

+ R v ' ' a i (t i )a i (t 2 } dt i dt 25 

m. (t + r 2 c 2 } j d2( . 

9 m U+RC -{{ ^ 


a l^ V dt l d V 


(4*36) 


We compare this with the variance of ( i b ( t)+i^( t)), given 
in eq. (4.22). 


Var [.ig) = 2k9RT 

Since R>> l/g m (normally, otherwise we would not have considered 
the inverse filter at all), the first term can be neglected. If 

)dt x dt 2 

(4.37) 

We can neglect this term. However if T is only of the order 
of the right hand side, we may not neglect this but approximate 
it as nearly 2K9RT. So now the net contribution due to 
thermal noise is 

Var [l g ] * 2K9RT x 2 » 4K9RT (4.38) 


T » 


RC 
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m 


T T 

IS 

o o 


d^ £ ( 

d Vt2 


V 


a l<V a l (t 2 
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This increase in thermal noise can be taken care of by 

modifying The decrease in the total thermal noise 

power due to high R is still present inspite of the doubling 

of Var[Z ]• We can see this more clearly by making an 
y 

analysis in frequency domain. 

The power spectral density of fi th (t) (due to i b (t) 

2k© 2 

and i A (t)) at the output the inverse filter S th * -jp • R 

= 2k© R 

The power spectral density of e^ ( t) (due 

to e A ( t) at the output the inverse filter S * «(R 2 c 2 w 2 +l) 

A OE y m 

where w = 2itf* 

These are as shown in Fig. 4.7. 



Fig. 4.7: p.s.d of © th (t) and e 0 (t) 

The total thermal noise p.s.d, is 

S T =S th+ S oE =2keR(l +5 i S+ ^ ■ (4.39) 

S T = S th + S oE = 2keR( 1+ RC 9 m 7 
If the bandwidth of the data signal is approximately l/T then 
the Bandwidth of the SS modulated signal is N/T. So for 
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the second term to be small conpared with 1 (i.e. if we 
want to be equal to N/T) 


RC 2 to 2 


I < 1 


ft 2 tcN 
oj = •' •■•■f 


i * © * 

RC 2 ( , ^i/T) 2 < 1 0 r T > VW5 m • 2 itNC (4.«) 

g ' m 

m 

if T = fR/g m .2nNC, then S o£ will be still less than S th but 
of the same order of magnitude as that of ^th* an 
alternative if we consider T to be such that both powers in 

S th and S w *° be ec * ual . Then, 


m N/T 9 9 

T = f 4 u> 2 df = |£ 2 . 4n 2 

0 9m ®m 


f 3 j N ^ T _ RC 


dE 2 /Nn3 
~ » 3 »v W 

y m 


or. 




RC£ 

g, 


4n d 


2tcNC 

TT- 


(4,41) 


- — i.e, T = VWi m • 

’m3 

resulting in a smaller value of T, 

To ensure that this equality of S th and S is maintained at 

the 0/P of the correlator also, we can choose the inverse 

filter characteristic such that it passes the signal frequency 

of interest only i.e. above f = N/T it attenuates the 

noises completely., 


Thus the receiver proposed solves the problem partially. 
The performance can again be found by assuming JL to be 
Gauss, r.v. Hence, 
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P e * erf c*( SNR 1 / 2 ) 

where the £NR now is slightly degraded than the first case. 

The variations of P g with optical power, <g>, R,N are sa»d as 
before (Chapter 5). The limitation of this receiver is in the 
last approximation we made regarding the noise contribution 
due to e^(t). For an amplifier* with low Sg only this is 
valid, otherwise not. Also arbitrarily large N is not 
possible due to the fact the p.s.d. of e^(t) increases as 
square of frequency and at higher values of frequency it tends 
to be' infinite, where upon the approximation will no more be 
valid. 

4.3 SIMULATION: 

The system considered in Section 4.1 is simulated on 
DEC-10* The details of the simulation are given in Appendix B. 
The results are presented in Chapter 5. The simulations- are 
done for no APD gain (i.e. <g>=l) since the probability 
density function of it is quite complex. 


CHAPTER 5 


PERFORMANCE RESULTS AND CONCLUSIONS 


This chapter begins with the presentation of the 
theoretical performance results for the receivers proposed 
in Chapter 4, The results are then interpreted and compared 
with results obtained by simulation . We then give a critical 
evaluation of the proposed system giving its merits and 
limitations. It concludes with some suggestions for further 
wo rk . 


5.1 NUMERICAL RESULTS ON THEORETICAL PERFORMANCE: 

As mentioned in Chapter 4, we take the performance 
criterion to be the probability of bit error given in terms 
of the signal to noise ratio SNR in eq. (4.27). 


Case I : Unlimited bandwidth availability: 

For convenience, we repeat here the expressions for 
the signal to noise ratio (SNR) and probability of error. 

P^ = erfc*(SNR X / 2 ) 


SNR = 


1.0 




+ <&) 2 . .^l/ft^/g ) 

11 e z <g> A bi ra 


v-2 , 
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Where all the symbols are as explained earlier. To 
get numerical results , we fix values fox some of these 
symbols. We choose to Oil assuming a silicon detector 
(for Germanium, K 0.5). & is chosen to be 300°K and the 

data symbol duration T is taken to be 1.562xl0~ 5 sec. This 
corresponds to data rate of 64K bps, typical of a digitized 
voice channel with 8 bit PCM. Use of ADM etc. can accomodate 
upto 2 voice channels in this data rate. Or it can even be a 
time multiplexed low data signals emanating from a user 
source of information. The amplifier trans conductance g m is 
taken to be 5m xj . The constants are 

hy/r) = S.XITxIO"' 3 - 9 !, k= 1.37xl0” 23 J/°K, e = 1.602xl0~ 19 coloumb. 

We have considered two cases of N=511 and 1023. The 
average value (worst case) C is dependent on the type of PN 
sequence and the code sequence length, N. Since M-sequences 
available for the above code lengths are only 24 and 30 
respectively. We did not consider them. However for 

comparison sake the multiple access interference term 

K 2 

Z a-,* fox K = 20 is tabulated in Table 5*1 for both 

mi 1,1 

M-sequences and Gold sequences. The Gold sequences axe 
considered for their large set of 513 and 1025 sequences for 
N=5ll and 1023 respectively. For these code sequence lengths 
the worst case is such that C/N is approximately 6'/.. We vary 
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th@ other parameters oyer a set values to observe the 

variation of the probability of bit error, P e « We first vary 

the photodetector (APD) gain and calculate the P . This is 

as e 

done for varying optical power as well£ varying K, the number 
of users. The variations are also observed with varying R, 
the equivalent bias resistance (from here onwards referred 
to as bias resistance ) and with both cases of N=511 and 
1023. The data of these are given in Tables 5. 2-5,7 and the 
plots are given in Figs, 5.1 to 5.6 where each figure corresponds 
to a specific N, and R. In all these cases we observe the 
performance first inproving then saturating, and finally 
deteriorating with increasing APD gain. 

We then evaluate the probability of error with varying 
optical power, again for certain values of K and for both 
values of N, The associated A?D gains are chosen to be 
optimum such that the p is minimum for: that value of optical 

V? 

power. These gains are obtained by solving eq. (4. 30) for each 
value of p. These variations are also observed with 
different R. The data are given in Tables 5.8 and 5,9 and 
the plots in Figs, 5.7 and 5.8. We then give the important 
variations of P G with K, the number of users in multiple 
access environment, again with optimum gains at chosen optical 
powers. This data is presented in Tables 5.10 and 5.11 
and plotted in Figs. 5.9 and 5.10. Vie finally present the 
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data in Table 5,12 and 5,13 of P versus R fox unit and 

c 

optimal APD gains. The plots are given in Figs, 5,11 and 
5.12. 

In the above case we assumed unlimited bandwidth at 
the receiver. But this is not so for high values of R due to 
non- zero capacitance values. For this problem the front end 
circuit is followed up by an inverse filter such that the 
signals are 'whitened’ as given in Chapter 4. 

The analysis of Section 2 of Chapter 4 indicates that 

in eq, (4.40 or 4.41) 

R,C,N, g m should be such that the condition^/ is satisfied for 
a given T to make the noise enhancement due to inverse 
filter small. However, T can be smaller than this resulting 
in a higher contribution from but still keeping the total 
thermal noise contribution less than what would have been 
obtained by keeping a low value of R in order to avoid the 
bandwidth limitations. As an Illustration we consider both 
the situations. 

1. R ts 1 -ft hence no ISI 

At the output of the amplifier total thermal noise spectral 
density (normalized by R^) = 2k©/R (1+ l/g m R^) 

S T » 201x2k9 for g m = 5 mX3~ 

2, R » 500 Kit and C = 10“ 11 F hence ISI is present. 
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At the output of the inverse filter, total thermal noise p.s.d. 

= . a' 2 + R- 2 . RW . 2S 

n 


g 


ym 


2kO (l/Rfl/g m R 2 + ) at f = N/T 

y m 


* 2k9x8 .4x10 


-4 


Thus introducing deliberate ISI by a high impedance 
front end, gives a considerably high performance over that 
obtained by low R front end circuit, even though there is 
enhancement of noise due to inverse filter. 


But for the above modif ica tions of the net thermal noise 
power, the system performance depends essentially the same 
way as it does in the first case. Due to the increase in 
this thermal noise power the performance slightly degrades at 
low optical powers. Thus it is appropriate to interprete the 
results of the previous case first and then see the changes 
required to get suitable performance in the present case. In 
the next section we try to interprete the results obtained with 
unlimited bandwidth case and compare them with the simulation 
results obtained. 
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5*2 INTERPRETATION OF THE PERFORMANCE RESULTS: 

First, we go through the variation of P e with APD 
gain. From Figs. 5.1 to 5.6, we see that 

1. The probability of error decreases first, then flattens, 

and finally increases, with increasing APD gain. This is 
so because of the fact (as given in eqn. 4,29) that the 
sum of shot noise and thermal noise decreases first with 
increasing APD gain. Thus the performance improves. But 
at medium gains the variance of the gain tends ^be such that 
the decrease in thermal noise is compen sated by an 
increase in the shot noise which tend to cancel out and 
performance remains constant. Further increase of the 
gain results in more shot noise contribution than thermal 
noise reduction leading to the increase of the total noise 
resulting in increase of P e . So there is an optimum value 
of <g> which minimizes the total noise which is found by 
solving * eq. (4.30). This result is consistent 

with the single user case. 

2. With increasing optical power the optimum gain decreases 
for a given number of users K. This is due to the reduced 
shot noise and thermal noise which is explicit in eq,(4.27). 
Thus a proper choice of the optical power has to be made 

to use as minimum APD gain as possible or even to cto 
a way with the APD. 
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3. With increasing K, the number of users, the optimum gain 
shifts left (or decreases). This is an important factor 
and is peculiar to the SSMA in fiber optic communication. 
This is due to the increase in number of users which results 
in a higher optical power (not that of signal) hence higher 
shot noise, unless the optical power is sufficiently high 
the shot noise power can exceed the multiple access power 
which normally exceeds the additive noise. 

4. Finally increasing R has an important effect on the APD 
gain. We observe the <g>^p^. decreases with increasing R. 

This is so because of the reduced thermal noise power due 
to increased R as we have seen earlier. Thus for FUlKSlthe 
optimum gain at p = 50 pW (for N=1023, K=50) is about 97 
where the corresponding value at ft=500Kfiis about 11 resulting 
in a performance difference of more than an order. Thus a 
proper choice of the bias resistance can obviate the need 
for a high gain APD. 

Then by observing the probability of error vs. optical power 
in Figs. 5.7 and 5.8 we can conclude the following. 

1. The probability of error is high at low optical power, 

falling off rapidly with increasing optical power. However 
such a decrease does not happen continuously. At high 
optical powers the effects of shot noise and thermal noise 
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be 

tend to^negligible and the multiple access interference term 

is dominant. Since this is independent of the optical power 

the probability of error saturates at a value which is 

‘irreducible*. This value can only be improved by choosing 

low cross correlation sequences and by taking a higher value 

of N, With doubling N, the bandwidth required also gets 

doubled. The use of an inverse filter is limited because 

(as shown in the previous section) the thermal noise power 

o 

enhancement is proportional to N . 

2. It is clear from these plots A an increase in R decreases 
the required optical power, to achieve a certain error 
rate It is more at low optical powers. The decrease 
is negligible at high powers due to the dominance of 
multiple access interference over thermal noise. But the 
low powers are only of interest to us where we can still 
have reasonable performance, with proper value of R . 
e.g. if the basic data signal with a rate of 64K bps 

-5 

is a voice channel, an acceptable error rate will be 10 . 

Thus at2^ iS value we see there is a decrease of the 

required optical power by around factor of four in going 

from a value of R of IKftto 500m. As mentioned earlier 

light 

these results are to be seen in £ of the peak optical 
power required. Thus at P e =10“ 5 we require the peak 
optical power to be (for K=50, N=1023) 2xkxp which is 
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equal to 40 nW at ffc=lK and 10 ny? at ft=500K. Since the 
saturation of the LED emitted optical powers can occur at 
high drive currents, thus making the system performance 
poor and making the relation between optical power and drive 
current nonlinear, a high value of R will always be preferred 
from this point of view. 

We then observe the variation of the probability with 
varying K at different optical powers with optimum APD 
gains. At low optical powers the dependence of P on K is 

V 

less due to the thermal noise and shot noise dominance. However 
still there will be an increase in P with increasing K 
because the shot noise power increases linearly with K, At 
high powers increases with increasing K due to increased 
multiple access interference. We have shown these cases for 
N = 511 and 1023 in Figures 5.9 and 5.10. It is then with 
increasing N the performance improves but at the expense of 
increased bandwidth have shown only up to K equal to 100 
because beyond this, the low performance obtained is not 
generally acceptable. 

We finally observe the variations of P e with R. The 
APD gain at these points aye again chosen to be optimal. At 
high optical power P Q is independent of R since the thermal 
noise power at both low and high values of R is negligible 
compared to the sicpal power. However at low powers, . ■ c.g. 
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at p = 1 nW for N = 1023 we obtain a performance at ft=500KA 

three orders of magnitude better than that obtained at 

^=10-0- ■ To get the performance of the former we need to 

increase the power by about 5 dB at R=10-A. Once again it is to be 

no ted, such an increase in optical power puts restrictions on 

the linearity of the LED. Thus the use of high impedance 

front end is justified at low optical powers. 

These conclusions will also be true for the limited 
bandwidth case but with slightly reduced performance. We also 
performed simulation (details in Appendix B) for the 
unlimited bandwidth case. Due to the complex nature of 
the APD gain probability density functions, we assumed unity 
gain (or PIN diode) case. Due to the heavy computation 
time required, we performed the simulation only for a low 
SNR case resulting in a theoretical P e of 2.5x10”^ with 
N=5ll, Using 10^ iterations the observed probability of error 
was about 1.5xlO"* 2 . Such a discrepancy could be due to the 
Gaussian approximation for the decision statistic and due 
to the independence assunption of shot noise on signal. From 
these facts some additional conclusions can be made about 
the receivers and the system as a whole. This is done in 
the next section. 
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5.3 CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK: 

Merits and Limitations of the Receiver:- 

The receivers in both the unlimited and limited bandwidth 
cases performed reasonably well at optical powers of practical 
feasibility. The particular problem of introduction of ISI is 
reasonably solved by the inverse filter. However, the net 
thermal noise reduction will only be smaller over that of the 
case where R is small when the code sequence lengtte axe 
of reasonable values. At very high code lengths the enhancement 
can be of the same order as or even more than that of the no 
ISI case* Then the use of high impedance front end is not of 
any use. The advantage we get at low values of code length 
is considerable from both the performance point of view as well 
as the less stringent requirement on the linear range of LED 
or Laser. As mentioned in [l] and observed here, in the 
asynchronous mode of operation, with SSMA, the number of new 
users with satisfactory performance is around N/lO. In a 
synchronous case this number will be much better. But such 
a restriction is what is not preferred and is in fact the 
major reason for choosing SSMA in preference to TDMA. So 
essentially SSMA proves superior to TDMA, by its asynchronous 
operation when the available bandwidth is more than the 
required bandwidth. In data networks, this is usually so 
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and SSMA can prove to be effective over TDMA. However the 
increased shot noise power and high peak power requirement 
of driving LED are some disadvantages of SSMA. 

The performance obtained by Gaussian approximation 
will be closer to the exact performance when 

i) the signal to noise ratio is relatively high 

ii) the number of asynchronous users is high which makes 
the multiple access interference tends to be Gaussian. 

Also the shot noise tends to be Gaussian when X(^) 
high most of the time [25]* This will be true for 
our high number of users case. 

iii) the APD gains are relatively low [23 ]* which is also 
true in our case. 

Suggestions for Further '/fork: 

In this thesis we considered only the suboptimal receivers 
for reasons of their simplicity and satisfactory performance. 
Unlike this, an optimum receiver structure may be derived to 
process the composite signal. In the system considered here 
we made an important assumption of the existance of perfect 
synchronization between the transmitted and locally generated 
PN sequences not only at the chip level but at the frame 
level also. In SSMA, such a synchronization is difficult to 
obtain and sometimes the number of users can be limited by 
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the synchronization problem only £24] , So it needs to be analysed 
how the performance changes with lack of perfect synchronization. 

The simulations are difficult to make for high data 
lengths due to the high CPU times required. A better way can 
be to inplement the system to check the validity of the 
performance results presented here. For this case synchronization 
can be directly made between transmitter and receiver. 

The case of unequal received powers of user signals, 
a case more realistic to SSMA is not considered here. However, 
such a situation can be considered which only modifies the 
multiple access interference. 

Finally, we also made the assumption of distortionless 
fiber channels. Though this might be true for the code 
sequence lengths we considered, at higher values of N the 
dispersion effects of the fiber can be analysed. Also the 
assumption of ideal nature of the LED or Laser (linearity and 
spectral coherence) and photo diode can be relaxed and the. 
non ideal case be analysed. 


Table 5.1 

Multiple Access Interference 
N = 511 K = 20 


, Interference value, 

of i x l0 3 /2 

M- Sequences 

Gold Sequences 

0.665 

0.641 

0.657 

0.640 

0.653 

0.718 

0.673 

0.662 

0.643 

0.716 

0.680 

0.678 

0.663 

0.726 

0.665 

0.607 

0.633 

0.671 

0.659 

0.662 

0.673 

0.684 

0.640 

0.659 

0.660 

0.694 

0.668 

0.644 

0.643 

0.653 

0.641 

0.652 

0.642 

0.681 

0.656 

0.676 

0690 

0.698 
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Table 5.2 

P e vs. Versus APD gain 
R = l.OKfl N = 511 


A?D gain 

P 3 

K = 10 

Probability of 
200 p W 

k « 20 

error, 

e 

p a 1000 pW 
K= 10 

K a 20 

1 

.44856E+00 

•44857E+00 

.25952E+00 

•26011E+00 

5 

. 25960 E+00 

. 26028 E+00 

• 11855E-O0 

.21088 E-0 2 

10 

.10076E+00 

•10374E+00 

.85783E-07 

.438 44 E-0 5 

15 

» 29664E-0 1 

. 33527 B-01 

, 31438 E-10 

•12223E-06 

20 

,69708 E-02 

.96947 E-02 

.16339E-12 

.19312E-07 

25 

•14088E-02 

.278 90 E-02 

.60600E-14 

.73248E-08 

30 

•26623E-03 

.87042E-03 

.77737E-15 

.43379E-08 

35 

.50914E-04 

,31128 E-0 3 

.21411E-15 

• 32678 E~08 

40 

.10521E-04 

.1C069E-03 

.95122E-16 

• 28264E-08 

45 

•24635E-05 

.64519E-04 

.57474E-16 

* 26580 E— 08 
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Table 5.3 

P Versus APD Gain 
ft = 91 KX1 , N « 511 


APD 

gain 


p==200pW 
K=10 K=20 


p = 500 pW 
K=10 K=20 


• 911E-01 .939E-01 .797E-03 .154B-02 
.734E-07 .571E-05 .336E-14. . 615E-08 
•114E-11 .896E-07 .177E-16 .177E-08 
. 327 E-l 3 .408E-07 .863E-17 .172E-08 
.130E-13 .410E-07 .947E-17 .202E-08 
•133E-13 • 522E-07 .131E-16 .245E-08 
.199E-13 .719E-07 .197E-16 .313E-08 
.347 E-l 3 .101E-06 .306E-16 .394E-08 
.647 E-l 3 .144E-06 .480 E-l 6 .495E-08 
•123E-12 .2Q3E-06 .754E-16 . 620E-08 


J 

p=100 pW 

K=10 K=20 

.3302-07 

•245E-05 

.304E-17 

• 756E*»09 

• 623E-18 

, 537E-09 

.593E-18 

• 576Er-09 

.712E-18 

.65XE-09 

.911E-18 

.744E-09 

.119E-17 

.853E-09 

• 157E-17 

.979E-09 

•207E-17 

.112E-08 

.273E-17 

•128E-08 


45 
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Table 5.4 

P Versus APD gain 
R = 50(3Ci2. N = 511 


Probability of error, (P 0 ) 

gain P~lG0pW p = 200 pW p - iooo pW 

K=10 K=10 K=20 K=10 K=2( 


1 

• 597E-01 

•144E-02 

5 

• 131E-07 

•641E-13 

10 

•462E-11 

•932B-15 

15 

.120E-11 

•585E-15 

20 

.163E-41 

• 167E-14 

25 

.350E-11 

•372E-14 

30 

•850E-11 

.849E-14 

35 

• 208E-10 

•191E-13 

40 

•493E-10 

.418E-13 

45 

.111E-09 

♦886E-13 


• 252E-02 

•242E-14 

.375E-08 

• 267B-07 

.389E-18 

.434E-09 

• 126E-07 

.357E-18 

.472E-09 

•165E-07 

.464E-18 

• 545E-09 

• 250E-07 

.621E-18 

•631E-C9 

. 38 5E-07 

.836E-18 

.7292-09 

• 588E-07 

.112E-17 

• 841E-09 

.883E-07 

.150E-17 

.969E-09 

.13PE-06 

.200E-17 

•111E-08 

.188E-06 

.266E-17 

•127E-08 
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Table 5,5 

P e Versus APD gain 



ft = IK _a 

N = 1023 

, K = 20 


APD Gain 

Probability of 

Error, P fi 



P = 50 pW 

P * 100 pW 

P = 500 pW 

P = 1000 

1 

• 487E+-00 

.474E+00 

.3736+00 

.299E+O0 

5 

.4356+00 

.373E+00 

.563E-01 

.124E-02 

10 

.373E+00 

• 26QE+0Q 

.132E-02 

.131E-06 

15 

• 314E+00 

.1696-00 

.162E-04 

.952E-16 

20 

.260E+00 

•103E+00 

.262E-06 

• 1016-11 

25 

• 2136+00 

.5996-01 

.900E-08 

.6966-13 

30 

• 171E+00 

.336H-01 

• 710E-09 

•149E-13 

35 

• 137E+00 

.1856-01 

.116E-09 

.6316-14 

40 

#1096+00 

» 102E—01 

.339E-10 

.402E-14 

45 

.8636-01 

.5826-02 

.1526-10 

•332E-14 
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Table 5.6 

P 0 Versus APD gain 
R as 91KJTL N = 1023 


APD Probability of Error, P 

Gain e 

P =» 100 pw p = 500 pW 

K = 20 K = 50 K = 20 K=50 


1 

• 250E+O0 

•251E+00 

• 84QE-03 

.211E-02 

5 

♦ 124&-02 

.403E-02 

.420E-13 

.315E-06 

10 

.155E-05 

.257E-03 

.911E-15 

.199E-06 

15 

.678 £-07 

• 1 47 E-0 3 

.833E-15 

.246E-06 

20 

.277E-07 

.163E-03 

.140E-14 

.332E-06 

25 

. 297E-07 

.219E-03 

.275E-14 

.452E-06 

30 

.467E-07 

.305E-03 

558E-14 

.612E-06 

35 

.844E-07 

.424E-03 

.112E-13 

. 8 20 E— 0 6 

40 

.157E-06 

.579E-03 

.222E-13 

*108E-05 

45 

291E-06 

. 773E-03 

.431E-13 

.142E-05 
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Table 5.7 

P G Versus APD gain 
R « 500KO. N = 1023 


APD 

Gain 


Probability of Error, p Q 
p « 100 pW p = 500 pW 

K=*20 K s 50 . K = 20 . K = 50 


1 

.603E-01 

• 661E-01 

.280E-08 

• 416E-05 

5 

• 12QE-06 

•842E-04 

.190E-15 

•117E-06 

10 

.114E-08 

.435B-04 

.197E-15 

•155E-06 

15 

•114E-08 

•652E-04 

.428 E- 15 

•221E-06 

20 

•271E-08 

• 107E-03 

.982E-15 

• 313E-06 

25 

.736&-08 

.171E-03 

•221E-14 

.436B-06 

30 

.192E-07 

.262E-03 

•482E-14 

• 596B-06 

35 

.467E-07 

.384E-03 

• 101E-13 

.807E-06 

40 

.104E-07 

•541E-03 

.206E-13 

.107E-05 

45 

•217E-06 

.736E-03 

.407E-13 

.140E-05 


% 

\ 
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Table 5.10 

P c Versus no. of users 


N 

* 511 R = 

500101 



K 

Probability of 

* 

Error, P e 



p = 100 pW 

p « 200 pW 

p = 1000 pW 

p=10000 pW 

10 

.138E-11 

• 996E-15 

.485E-18 

.927E-19 

20 

.353E-06 

.132E-07 

.463E-09 

•223E-09 

30 

* 209E— 04 

.257E-05 

.320E-06 

.203E-O6 

40 

.185E-03 

.397E-04 

.867E-05 

. 623E-05 

50 

•700E-03 

.206E-03 

. 622E-04 

.479E-04 
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Table 5.11 
Versus no. of users 

N = 1023 R * 500K jQ. 


Probability of Error, P 

e 

p=500 pW p = 1000 pW p = 10,000 pW 


20 

.155E-15 

.686E-17 

.409E-18 

30 

.126E-10 

.183E-11 

• 322E-12 

40 

.405E-08 

.99 E-09 

.280 E-09 

50 

.120E-06 

.396&-07 

• 145E— 07 

60 

.118E-05 

.480E-06 

.213E-06 

70 

•625E-05 

•291E-05 

♦145E-05 

80 

#21 6E-04 

•111E-04 

.609E-05 

90 

•563E-04 

.313E-04 

f 184E-04 

100 

•124E-03 

. 735E-04 

• 456E—04 



Table 5.12 : p versus bias resistance 
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Table 5.13 

P 0 Versus Bias Resistance 


K. ** 

R, 

XL 

1023 

K. 

= 50 


Probability of Error, P g 
p=100 pW p=1000 pW p = 100 pW p= 1000 pW 

10 

0 • 499&f00 

0*493&f00 

0.208E-00 

0.590E-07 

100 

0.494E+00 

0.448E+-00 

0.585E-01 

Q.208E-07 

1000 

0.474Bf00 

0.260E+00 

0.809E-02 

0.163E-07 

10000 

0.412E+00 

0.198E-01 

0.992E-03 

0.152E-07 

100000 

0.241E+00 

0.742E-05 

0.135E-03 

0.147E-07 

500000 

0.661E-01 

0.114E-06 

0.426E-04 

0.145E-07 

1000000 

0.210E-01 

0.522E-07 

0.278E-04 

0.145E-O7 
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APPENDIX A 


SELECTION AND GENERATION OF GOLD SEQUENCES 

As we have seen in Chapter 3, Gold sequences are large 
in number for a given code sequence length, N and also 
possess low cross correlation properties. For a given N> they 
can be selected as follows. 

Let f^x) be the minimal polynomial of the primitive 
element oc€GF(2 n ) and f^.(x) be the minimal polynomial of 
where f^(x) and f ^( x) are of degree n and if n and t are 
such that 


= 2 n+1 ^ 2 + 1 , 

n odd 

= 2 n+2/2 + 1, 

n even 


then the product f ( x) == f^( x.)f ^.( x) determines a linear 
feedback shift register (LFSR) which generates the 2 n +l 
different Gold sequences (corresponding to the 2^+1 states 
in distinct cycles) of period (or length) N = 2 n -l[l], To 
state this in terms of the A4- sequences generated by LFSRs 
corresponding to f^(x) and f^(x) let u be the M-sequence 
generated by an LFSR corresponding to f^(x) and v be the 
similar JA- sequence Corresponding to f^(x). Then the set of 
Gold sequences G(u r v) is given by [22] 

G(u,v)r|u,v, u © v, u0Tv, u (P) T^u © T 3 v, . „ . u©T^ 
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where 0 indicates term wise modulo-2 addition of the two 
sequences and rSa indicates the ith cyclic shift of the 
sequence u>. It is to be noted that in the above set, except 
for u and v any other sequence is not an M-sequence, i.e. 
they do not have the properties of JvU sequences. The generation 
of Gold sequences is clear from the above. This can be 
generated in 2 ways. One is to have an LFSR corresponding 
to f ( x) and generate the different sequences by inserting 
different initial states into the LFSR. The other method is 
to have two LFSR 1 s corresponding to f^x) and f^x) and take 
the modulo 2 sum of their outputs. Wo have indicated both 
the methods in Fig. A.l for 

f^Cx) = l+x 2 +x 5 and f^Cx) = fg(x) = l+x 2 +x 4 +x 5 , 
since t s _ g 

The polynomials are commonly represented, in octal notation. 

Then f ^ C x) is represented as 25(=sl0010l) and f^(x) as 
65(=1!0101) . 

In our case we choose [26] the polynomials for n=9,10 

i 

corresponding to N=511 and 1023 respectively as given 
below in octal notation. 

N n t f x (x) f t (x) 

511 9 33 1021 1333 

65 2011 3515 


1023 


10 
















APPENDIX B 


COMPUTER SIMULATION 

As stated in Chapter 4, we made simulations for the 
unlimited bandwidth receiver. The simulations are made for 
only unit APD gain because of the difficulty in generating 
the random gains vdth its probability density functions which 
is quite complex. We know the intensity A( t) is given by 

qb c K 

A(t) a — (K + S d.(t~ ^,)a (1> T.)) ( B.l) 

hy T i=l 1 11 1 

The term inside the brackets changes between O and 2K. 
Also It changes, say for t£L(0,T), possibly for each of the 
chip intervals. So to simplify the situation we assumed 'Y^ 
to be taking only discrete values which are multiples of T c , 
the chip duration. This makes ,\(t) to have a fixed value 
during each .chip interval, through randoms because of the 
randomness of 2<i<K. We removed the i=l case because we 
assume the data bits of user 1 signal to be all ones. Such 
an assumption cbes not lead to any inferior performance. The 
simulation is performed as follows: 

1. The random delays s are generated using uniform distri- 
bution for integers between zero and N. 

2. The data bits d. ,, d. are generated for each i assuming 

1 f 1 f o 

the correlation interval to be in (o,T). The PN- sequences 
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a j[(t), 2 < i <}C arc modified by these 'IT’S and d^’ s and 
.added to a^(t) to get the ^X(t) for t£(o,T). 

3. Using thiS/\(t) we generate N Poisson random variables 
corresponding to each chip interval. These Poisson random 
variables represent the total number of electrons emitted 
in that chip interval. The emission time instants are not 
important to us because the correlating PN sequence has 
constant value (or it integrates over the chip interval) in 
each chip interval. So these Poisson random variables 
are weighted by the corresponding N chip values of a^(t) 
and 1 summed up to get the shot noise signal contribution 
to the decision statistic. 

4. We then generate a Gaussian random, variable and add this 

to the above quantity to make the total decision statistic. 
We compare this statistic to the threshold and decide as 
the data bits of user 1 signal. 

Since detection of each data bit needs the generation of N 
Poisson random variables, the computation time required is very 
high. For a data length of 1000 bits the confutation time 
required is about one hour for N=511. So higher values of 
N, are not considered. 

The program listing of the simulation is 
given after the list of References. 
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